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Abstract 

Quantum discord is a general measure of bipartite quantum correlations with a potential role in 
quantum information processing tasks. Spin clusters serve as ideal candidates for the implementation 
of some of the associated protocols. In this paper, we consider a symmetric spin trimer and a 
tetramer which describe a number of known molecular magnets and compute the quantum discord 
in the ground and thermal states of the clusters. The variations of the quantum discord as a 
function of an anisotropy parameter, magnetic field and temperature are investigated. We obtain 
a number of interesting results such as a finite value of the quantum discord in the trimer ground 
state for which the pairwise entanglement is known to be zero, differences in the nature of some 
of the variations in the ferromagnetic and antiferromagnetic cases and discontinuous jumps in the 
magnitude of the quantum discord at first order quantum phase transition points. A remarkable 
feature that is observed is that the quantum discord completely vanishes only in the asymptotic 
limit of temperature T — > oo. We further study the dynamics of the quantum discord and the 
pairwise entanglement at T = under the effect of a dephasing channel describing the interaction 
of the reduced spin cluster state with independent local environments. The QD is found to vanish 
asymptotically as t — > oo. In the case of the spin trimer, the pairwise entanglement has a zero value 
at all times and reaches a zero value in a finite time in the case of the tetramer. 



PACS Numbers: 03.67.-a,03.65.Ud,03.65.Yz,75.10.Jm,64.70.Tg 



1 Introduction 

An interacting quantum system is characterized by the presence of correlations amongst its different 
constituents. The correlations have, in general, both classical and quantum components. The most 
prominent example of quantum correlations is that of entanglement which serves as the fundamental 
resource in several quantum information processing tasks such as quantum computation, teleportation 
and dense coding [lj. Entanglement can be of various types, e.g., bipartite, multipartite, zero temper- 
ature, finite temperature etc. for which a number of quantitative measures exist. The entanglement 
content of quantum states and its variation as a function of changing parameters have been extensively 
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investigated in recent times [2,[3j|I]. A different measure of pairwise quantum correlations, namely, the 
quantum discord (QD) has further been proposed based on the information-theoretic concept of mutual 
information j5j El [7J El |9]. The basic difference of the QD from entanglement is evident from the fact 
that the QD is non-zero in certain separable states which are, by definition, unentangled. The QD is 
defined to be the difference between two quantum extensions of the classical mutual information. In the 
classical domain, the two representations are exactly equivalent. 

In classical information theory, the correlation between two random variables A and B is measured 
by their mutual information [I] 

I(A,B) = H(A)+H(B)-H(A,B) (1) 

The random variable A takes on the values 'a' with probabilities given by the set {p a }- H(A) = 
— ^ a p a log 2 p a is the Shannon entropy. H(A,B) corresponds to the joint Shannon entropy defined as 
H(A, B) = — YlabPa-yb ^°S2Pa,b- As alternative representation of the mutual information is given by [HI [9] 

J (A, B) = H(A) - H{A\B) (2) 

where H(A\B) is the conditional entropy and quantifies the lack of knowledge of the value of A when that 
of B is known. The exact equivalence of the expressions in equations ([1]) and (J2J) can be demonstrated 
using the Bayes' rule pQ. 

The generalisation of the expressions to the quantum case is achieved via the replacement of the 
classical probability distribution and the Shannon entropy by the density matrix p and the Von Neumann 
entropy, S(p) = —tr (plog 2 p), respectively. Thus, the quantum versions of equations ([T]) and (j2J) can be 
written as 

I (pab) = S {pa) + S (p B ) - S (p AB ) (3) 
J (pab) = S (p A ) - S (p A \pb ) (4) 

where S (pab) is the quantum joint entropy and S (pa\Pb) the quantum conditional entropy. The 
latter quantity is, however, not properly defined by a simple replacement of the Shannon entropy by 
the Von Neumann entropy. The magnitude of the quantum conditional entropy, by the very nature of 
its definition (ignorance of A once B is known), depends on the type of measurement. Since different 
measurement choices yield different results, equations ([3]) and are no longer identical. We consider 
Von Neumann-type measurements on B defined in terms of a complete set of orthogonal projectors 
{Hi}, corresponding to the set of possible outcomes i. The state of the system, once the measurement 
is made, is given by 

Pi = (/®nf)PAfl(/®nf)/ Pi (5) 

with 

p i = *r((/®nf)p AB (/®nf)) (6) 

/ denotes the identity operator for the subsystem A and pi gives the probability of obtaining the outcome 
i. From equation HI an alternative expression of quantum mutual information is given by [8, 9j 

j {pab, {nf }) = s( PA ) - s ( PAB |{nf }) (7) 

The quantum analog of the conditional entropy is 

S{p AB \{Uf})=^2p i S{p i ) (8) 



2 



Henderson and Vedral (7] have shown that the maximum of J (pABi {nf }) w.r.t. {11,,} provides a 
measure of the classical correlations, C (pab), i-e., 

C{p AB ) = max(J{p AB ,{Bf})) (9) 
{nf} V 

The difference between the total correlations I (pab) (equation [3]) and the classical correlations C (pab) 
defines the QD, Q (pab), 

Q (p AB ) = I (Pab) -C ( P ab) (10) 

In the case of pure states, one can show that the QD reduces to the entropy of entanglement [9] so 
that entanglement provides the sole contribution to quantum correlations. In the case of mixed states, 
however, the QD and entanglement provide different measures of quantum correlations. QD has been 
quantified in a number of two qubit states (HI El HQj [HJ [12]. For spin Hamiltonians with certain 
symmetries, the two-spin reduced density matrix p^ (the two spins are located at the sites i and j) can 
be expressed in the basis {| ff), | fj,), | ||), | H)} as 



Pij 



( a / \ 

o h z o 

z b 2 

V / d J 



[ID 



The elements of the reduced density matrix can be expressed in terms of the single-site magnetization 
and two-spin correlation functions. As shown in Ref. [8], the eigenvalues of pij are 



A = 7<! (l + c 3 ) + a/(c4 + c 5 ) 2 + (ci-c 2 ) 2 



Ai = ^ (1 + c 3 ) - J (c 4 + c 5 ) 2 + (ci - c 2 ) 2 



A 2 = 7 <{ (1 - c 3 ) + a/ (c 4 - c 5 ) 2 + (ci + c 2 ) 2 



A 3 = -<j(l-c 3 )- \/(c i -c 5 ) 2 + (c 1 + c 2 ) 2[ > (12) 



with 



ci = 2z + 2/ 

c 2 = 2z-2f 

c 3 = a + d — bi — b 2 

C4 = a — d — bi + b 2 

c 5 = a — d + bi — b 2 (13) 
The mutual information (equation [3J can be written as [SJ [H] 

3 

I (pab) = S (p A ) + S (p B ) + J2 X « 1o S2 A « (14) 

o=0 



3 



where 



o / x (1 + c 5 ) . (l + c 5 ) (1 - eg) , (1 - c 5 ) 
5 (pa) = = log 2 = = log 2 z — 
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C / \ (I + C4), (I + C4) (I-C4) (I-C4) , , 

S(Pb) = log 2 log 2 (15) 

The reduced density matrix pij (equation [TTj) has a simpler form when specific symmetries of the spin 
Hamiltonian are taken into account. The element / = when the ^-component of the total spin com- 
mutes with the Hamiltonian, i.e., is a conserved quantity. Also, a — d, b\ — 62 when the magnetization 
density has expectation value zero resulting in C4 = C5 = and C\ = c 2 (equation [T5|) . Under these 
simplifications, the maximization procedure for computing the classical correlations C (pab) (equation 
[9]) can be carried out analytically to yield [8j [9] 

C (Pab) = ^y^ log 2 (l - c) + log 2 (l + c) (16) 

where c = max(|ci| , |c 2 | , |c 3 |). The QD, Q (pab), is given by [8, 9J 

Q {pab) = I {pab) - C (pab) 

= -[(1 - Ci - 02 - c 3 ) log 2 (l - Ci - c 2 - C 3 ) + (1 - Ci + c 2 + c 3 ) log 2 (l - Ci + c 2 + c 3 ) 
+ (1 + Ci - c 2 + c 3 ) log 2 (l + Ci - C 2 + C 3 ) + (1 + Ci + c 2 - c 3 ) log 2 (l + Ci + c 2 - c 3 )] 

log 2 (l - c) - ^1 i og2 (i + C ) (17) 



The ground and thermal state entanglement properties of small spin clusters have been computed 
in earlier studies [T5J HU H5]- in this paper, we compute the mutual information, classical correlation 
and QD in the ground and thermal states of a symmetric spin trimer and a spin tetramer with nearest- 
neighbour (n.n.) as well as next-nearest-neighbour (n.n.n.) interactions. There are several examples of 
molecular magnets represented by spin trimers and tetramers j 16] . In sections II and III, the results 
for the spin trimer and the tetramer respectively are presented. In section IV, we analyze some earlier 
results p2] on the robustness of QD to sudden death in the context of the spin trimer and tetramer. 
Section V contains a summary of the main results obtained in this paper and concluding remarks. 



2 Classical and Quantum Correlations in Spin Trimer 

The symmetric spin trimer consisting of three spins of magnitude | is described by the Heisenberg 
exchange interaction Hamiltonian 



4 



where Sf (a = x, y, z) defines the spin operator at the ith site of the trimer, e is an anisotropy parameter 
(e < 1) and J is the strength of the exchange interaction. The eigenstates of the trimer are given by 



1) 

2} 

3) 

4) 

5) 

6) 

7) 
8} 



I ttt) 

^(gim> + g 2 i nt) + utt» 
^(<? 2 im> + <?i t;t> + utt>) 
^= (I m) + 1 tit) + utt» 



(q\nt) + q 2 \m) + \m)) 



-^{q 2 \m) + q\m) + \m)) 



s/3 



(i lit) + 1 m) + 1 m» 



(19) 



where q = e 1 ^ is the cube root of unity satisfying g 3 = 1 and q + q 2 + 1 = 0. The corresponding energy 
eigenvalues are 







F 3J 




E 2 




(l + 2e)J 

E 3 = E Fj = E 6 = 




E 4 




(l-4e)J 
^= 4 


(20) 


We first calculate the ground state 


(T 


= 0) QD in the isotropic case e = 1. 


There are two distinct 


eigenenergies given by 








ei 




E 1 = E 4 = E 7 = E 8 = — 




e 2 




E2 = E 3 = E 5 = E 6 = — — 


(21) 



In the case of antiferromagnetic (AFM) exchange interaction, J > 0, the ground state is four-fold 
degenerate with the energy e%- The ground state density matrix is given by 



AFM 



(|2)(2| + |3)(3| + |5)(5| + |6)(6|) 



r 9 4 

The two-qubit reduced density matrix in the standard basis has the form shown in equation [TT1 with 



(22) 



a = d = -, b\ 
6' 

From equation [131 c = max (l c i| > l c 2| , | c 3 1 ) 
(equation [T2|) with c 4 = 0, c 5 = 0): 



D 



(23) 



|. The eigenvalues of the reduced density matrix are 



Ai — A2 — A3 — — , A4 
o 



1 
2 



(24) 



From equations (JHJ), (1T51) . ( TIE]) and ( 1TT1) . the mutual information / (pab), the classical correlation 
C (pab) and the QD, Q (pab), are: 

/ (pab) = 0.207, C (pab) = 0.082, Q (p AB ) = 0.125 (25) 

We next consider the ferromagnetic (FM) case with J < 0. The ground state is four-fold degenerate, as 
in the AFM case, with the ground state energy e\ (equation I2TT) . The ground state density matrix p^ M 
is given by 

^ = ^(|1><1| + |4>(4| + |7><7| + |8><8|) (26) 
The reduced density matrix has the same form as in the AFM case with 

a = d = \,b l =b 2 = \, z = \, f = U (27) 
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The value of c is | and the eigenvalues of the reduced density matrix are 

Ai = A 2 = A 3 = ~ A 4 = (28) 

The values of the mutual information, classical correlation and QD are 

/ (p AB ) = 0.415, C (pab) = 0.082, Q (p AB ) = 0.333 (29) 

From equations ( 1251) and ( |29|) . one finds that the magnitude of the classical correlation C (pab) is the 
same in the FM and AFM cases whereas the QD has a lower value in the latter case. 
We next compute the thermal state QD. From equation [201 the partition function is 

Z = 2 (e-^ 1 + 2e (1+2eh + e (1 " 4eh ) (30) 

where 7 = ., J r . The thermal state density matrix is 



4fc s T 

Z 



l>r =iEe-^|i><<| (31) 



i=l 

with (3 = -j^-f- The reduced thermal state density matrix has the form shown in equation [TT1 with 



bl = fe 2 = A ( 2e d+2eh + e (i-4 £)7) 



z 



2i_ ( e (l-4e) 7 _ e (l+2e) 7 ^ 



3Z 

f = (32) 

We calculate C (pab) and QD as functions of the parameter e and temperature T. Figure [T] shows 
the variation of the QD and C (pab) (inset) with the anisotropy parameter e. In both the AFM and 
FM cases, the QD increases with e. At a fixed value of e, the QD decreases with increasing T in the 
AFM case though non-monotonic behaviour is observed in the FM case. Furthermore, C (pab) (inset) 
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0.14 




Figure 1: Variation of quantum discord (QD) and classical correlation C (pab) (inset) as a function of 
the anisotropy parameter e for different temperatures in the (a) AFM and (b) FM cases with \J\ = 1. 




0.0 0.2 0.4 0.6 0.8 1.0 



Figure 2: Variation of quantum discord (QD) and classical correlation C (pab) (inset) as a function of 
temperature T for different values of the anisotropy parameter e in the (a) AFM and (b) FM cases with 
|J| = 1. 
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increases slowly with e for a fixed value of T in the AFM case whereas in the FM case it decreases with 
e. In both the cases, C (pab) decreases with T for a fixed value of e. 

Figure [2] shows the variation of QD and C (pab) with temperature for different values of e in both 
the AFM and FM cases. In the AFM case, the QD decreases with T for both low and high values of e. 
In the FM case, for low values of e, the QD first increases with T, reaches a maximum value and then 
decreases. For high values of e, the same features as in the AFM case are observed. For the AFM trimer, 
it is well-known [18] that there is no pairwise entanglement both at T = and at finite temperatures. 
For the FM trimer, the same is true when < e < 1. The calculations of the present study show 
that the QD, a different measure of quantum correlations, has non-zero values both in the ground and 
thermal states. Another remarkable feature is that though the magnitude of the QD decreases as T 
increases, it falls to zero only asymptotically, i.e., as T — > oo. This is in contrast to the property of 
pairwise entanglement in the thermal states of spin clusters |T4"1 H~5l ITS] . The concurrence, a measure 
of the pairwise entanglement, falls to zero value above a finite temperature. A proof that the QD has 
a non-zero value in the thermal states of the symmetric trimer over an extended temperature region is 
obtained using the result quoted in Ref. [19] . The two-qubit reduced density matrix (equation [TT]) has 
vanishing discord if the following conditions are satisfied: 

\f \ = \z\, a = b 2 , d = bi (33) 

The reduced density matrix elements for the spin trimer are given in equation [32] and one can verify 
that the conditions set in equation [33] are obeyed only when T — > oo. 

We now briefly consider the case when an external magnetic field term, h^^ i=l S*, is introduced in 
the spin trimer Hamiltonian (equation [IB]) . The energy eigenvalues are now 

3J 3h 

E x = 1 

1 4 2 



I 

E A = 



h (l + 2e)J 

2 4 
h (1 - 4e) J 



E 5 = E t 
E 7 = 



2 4 

h (1 + 2e) J 



2 4 
h (l-4e)J 



2 4 

= ^-f (34, 

One has a ^ d in equation [TT1 so that C4 = C5 7^ (equation [TBI. In this case, QD has to be computed 
numerically. We briefly discuss the results obtained for the case of e = 1. In the AFM case, the ground 
state is doubly degenerate for h < ^ and the ground state density matrix is 

^=~(|5)<5| + |6}<6|) (35) 
The elements of the reduced density matrix (equation [11]) are 

a = 0, d = ~, h = b 2 = \ ) z = -\ ( 36 ) 
3 3 o 

3 J 



The numerically computed value of the QD is 0.125815. For all values of h > 4jr, the ground state is the 

3J 
2 



fully separable state |8) with QD = 0. At h = a first order quantum phase transition (QPT) takes 



0.15 



0.10 

Q 

a 

0.05 

0.00 I mmh 

0.0 0.5 1.0 1.5 2.0 

h/J 

Figure 3: Variation of quantum discord (QD) with magnetic field h at T — 0. The QD remains constant 
in the range < h < At h = ^, the value of QD is 0.0838764. When h is > the QD is zero as 
state 1 8) is the ground state. 
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Figure 4: Variation of quantum discord (QD) as a function of temperature for different values of h and 
e = 1.0 (| J\ = 1.0) in the AFM and FM (inset) cases. 

place and the ground state is three fold degenerate: 

P 9 ,' QPT = ^(|5)(5| + |6)(6| + |8)(8|) (37) 

The reduced density matrix has the elements 

5 2 1 

o = 0, d = g> h = b 2 = -, z = -- (38) 

The QD has the value 0.0838764. Figure [3] shows the variation of QD with the magnetic field h at 
T = 0. Figure H] shows the variation of QD with h in the AFM and the FM cases. For comparison, the 
plots for h = are also included. In the AFM case, the effect of h on QD is not prominent. In the 
FM case, the QD decreases as h increases. For h — 1, the QD is zero at T = 0, rises with temperature 
initially and then assymptotically goes to zero. 

3 Classical and Quantum Correlations in Spin Tetramer 

The spin tetramer cluster is shown in figure [5j The n.n. spins interact via the isotropic Heisenberg 
exchange interaction Hamiltonian with interaction strength J\. The n.n.n. spins interact with exchange 
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interaction strength J 2 . We consider the AFM case with Ji, J 2 > 0. 



The eigenstates of the Hamiltonian are given by 

ii) = | trrt) 

|2> = ^(itm)+*imt)Htm)-mttt)) 
|3> = -^=(ittt4>-<imt>-mtt>+iuttt» 
|4> = (i tm> + 1 mt) + 1 tut) + 1 mt)) 
1 5} = 4= (i tm> - 1 mt> + 1 tut) - 1 tttt» 
16) = -L(|ttu)-iutt)+iimt)-ium)) 
|7) = -±=(\nn)-\nn)-i\mv+i\mi)) 
|s> = ^= (I tm) - 1 mt)) 

|9) = ^(ittm + ltm) + imt) + imt) + itm) + imt)) 
I io) = -j= (i tm) + 1 wtt) - 1 tut) - 1 mt)) 
in) = -1=(2| tm> + 2 1 mt) - 1 mt) - 1 mt) - 1 mt) - 1 mi)) 

V 12 

|i2) = ^=(imt)+mm)-imi)-^i tm)) 

1 13) = -J=(imt)-mm)-imi)+^i tm)) 

|i4) = ^(imt) + um) + imt) + itm)) 

lis) = ^(imt)-imt) + imt)-itm)) 
1 16) = imt) 

10 



The corresponding eigenvalues are 



ei 


= E 1 = 


Ei = Eg = Eu = 


= #16 = 




Co 


= £o = 


C/q = frjCL — /^v7 — 




Q — 


e 3 


= £5 = 


#8 = #15 = — ^1 


2 




e 4 


= i? 10 = 


3J 2 
2 






e 5 


= E n = 


= -2J, + # 







2 

2 



(41) 

We first consider the case J% > J 2 . The ground state is given by |11), which is an example of a resonating 
valance bond (RVB) state [20j HH |2T| . The RVB state is a linear superposition of two valance bond 
states in one of which there is a pair of horizontal valance bonds whereas in the other the valance bonds 
are vertical. A valance bond represents the singlet spin configuration ^ (| yj,) — I it))- The elements 
of the n.n. and n.n.n. reduced density matrices are 

1 ^ 1 

nn inn l.nn inn nn enn n 

a = d =-A =b 2 =-,z =--,/ =0 

a nnn = jnnn = 1 yann = = _ /n n = 1 jnnn = q /^n 

3 6 6 

The corresponding classical and quantum correlations are 

C nn ( P ab) = 0.350, Q nn (p AB ) = 0.442 

C nnn ( P ab) = 0.082, Q nnn {pab) = 0.333 (43) 

When Ji is < J 2 , the ground state is given by the state 1 10) which is again a RVB state |14| . The n.n. 
and the n.n.n. reduced density matrices have the following elements 

a nn = jnn = _ fe nn = fo nn = _ ? /n __ pn _ g 
a nnn = jnnn = g^ yvan = yvan = _ ^nn = _^ ynnn = g (44) 

2 2 

The corresponding classical and quantum correlations are 

C nn (PAB) = Qnn (PAb) = 0.0 
C nnn {pAB) = Qnnn (PAb) = 1-0 (45) 

At Ji = J 2 , a first order QPT takes place with the ground state density matrix given by 

p^=~(|10)(10| + |ll)(ll|) (46) 

The elements of the n.n. and n.n.n. reduced density matrices are 

1 11 

a = d = -, b 1 = b 2 = -, z = -- (47) 
6 3 6 
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Figure 5: A tetramer of spins of magnitude |. J\ and J 2 denote the strengths of the n.n. and n.n.n. 
(diagonal) exchange interactions. 

The classical and quantum correlations are 

C nn (pab) — C nnn (pab) = 0.082 

Qnn(pAB) = Qnnn (Pab) = 0.125 (48) 

We next calculate the classical and quantum correlations in the thermal state of the tetramer. The 
partition function of the system is given by 

Z = 5e~ /3ei + Qe~ Pe2 + 3e~^ 3 + e~^ 4 + e"^ 5 (49) 

where the energies e^'s are as listed in equation HU The n.n. reduced density matrix has the elements 

a nn = d nn = }_ -fo + 3 _^ + 1 + I ^ + l_ fi 

/j \ o Z Z 4 LZ 

b nn = h nn = }_( 5 -f)e x + 3 -/»ea ^ 3 1 ^ 1 -/3e 5 

1 2 Z V6 2 4 12 



Z \6 2 3 

f nn = (50) 

The n.n.n. reduced density matrix has the elements 

a nnn = jnnn = 1 / + e -/3e 2 + e ~l3e, + 

Z V 3 3 

6?™ = &o nn = - ( -e"^ 1 + 2e- /3e2 + ie"^ 3 + -e"^ 4 + -e~ pe ' 
1 2 Z \6 2 2 6 

1 / 5 _ft„„ 1 __«o„ 1 1 



Z \6 2 2 6 

f nn = (51) 

Figure [6] shows the variation of n.n. concurrence (CN), QD and classical correlation (CC) with temper- 
ature for the AFM case. When J x > J 2 , both the CN and QD decrease with temperature but QD has 



12 




Figure 6: Variation of quantum discord (QD), concurrence (CN) and classical correlation (CC) as 
functions of temperature for (a) J\ > J2 {J\ = 2J 2 = 1.0) and (b) J\ < J 2 (2J\ = J2 = 1.0). The results 
correspond to the n.n. reduced density matrix. 
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Figure 7: Variation of quantum discord (QD), concurrence (CN) and classical correlation (CC) as 
functions of temperature for (a) J\ < J2 (2Ji = J2 = 1.0) and (b) J\ > J2 {J\ = U2 = 1.0). The results 
correspond to the n.n.n. reduced density matrix. 

non-zero values at temperatures much higher than the value at which CN becomes zero. When J\ < J2, 
CN has zero value at all temperatures. The QD is zero at T = 0, then it increases with temperature 
to reach a maximum value after which it decreases with temperature. Figure [7] shows the variation of 
the n.n.n. CN, QD and CC with temperature in the AFM case. The CN has non-zero values only 
when Ji < J 2 . The magnitude of CN decreases with temperature and falls to zero value at a specific 
temperature. As in the n.n. case, the QD has non-zero values at much higher temperatures. 

4 Quantum Correlations under Decoherence 

The interaction of a quantum system with its environment results in decoherence, i.e., a destruction 
of the quantum properties including correlations of the system. The dynamics of the two-qubit QD 
under different types of environment have recently been investigated [T7] |22| [23| |24"] . In the case of 
local environments, each qubit interacts with its individual environment. The channel representing the 
interaction between a qubit and its environment can be of various types: amplitude damping, dephasing, 
bit flip, phase flip etc. [Ij. In this section, we study the dynamics of the two-qubit entanglement and 
QD under the influence of a dephasing channel. The initial (time t = 0) two-qubit state is the Werner 
state 

p(0) = (l-a)- + ahrWI (52) 

with a G [0, 1], / the identity matrix and \ip~) = -h^ (| ti) — I it))- We note that the two-qubit reduced 
density matrices of the ground states of the AFM spin trimer and tetramer represent the Werner states 
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Figure 8: Variation of quantum discord (QD) and concurrence (CN) with 7 for a = | (equations ([5 
and (158)) ). The two-qubit reduced density matrix is obtained from the AFM ground state of the trimer 



with a = I (trimer) and a = | (tetramer, J\ > J 2 ) respectively. The dynamics of a two-qubit state 
under the effect of the dephasing channel and with the Werner state as the initial state have been studied 
in Ref.|17|. We utilize the results of this study (some minor errors have been corrected) to examine the 
evolution of the QD in the two-qubit states with a = | and |. 

In the Kraus operator representation, an initial state, p(0), of the qubits evolves as [T7] 

p{t) = Y J E,Mml,u (53) 

where the Kraus operators E^ v = E^® E v satisfy the completeness relation ^ v E^EJ^ v = I for all t. 
In the case of the dephasing channel, the Kraus operators have the matrix form 

where 7 = 1 — e~ Tt with V denoting the decay rate. The elements of the density matrix of the two-qubit 
system, with p(0) given by equation [52] evolve to 

Pu{t) = Pu(0), i = 1,...,4 (55) 

P23(t)=P2 3 (0)(l-7)=P32(t) (56) 

The concurrence for the two-qubit evolved state is 



cxin.w)-- (^-7)-^ (57) 



The QD for the state is given by 



Q (pab) = \ {F(a + b) + F(a-b)}-^P- (58) 



with F(x) = xlog 2 x, a = (1 + a) and b = 2a(l — 7). Figure [8] ([9]) shows the variation of CN (pab) 
and Q (pab) as a function of 7 for a = | (|). In the case of the tetramer, the concurrence becomes zero 
in a finite time whereas the QD vanishes only in the asymptotic limit t — > 00. Thus, the QD, unlike 
quantum entanglement, exhibits robustness to sudden death. 
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Figure 9: Variation of quantum discord (QD) and concurrence (CN) with 7 for a = | (equations (157|) and 
(|58|) ). The two-qubit reduced density matrix is obtained from the AFM ground state of the tetramer. 

5 Discussion 

In this paper, we have calculated the QD of the two-qubit ground and thermal states of a symmetric 
spin trimer and a tetramer with both n.n. and n.n.n. exchange interactions. In both the cases, the 
QD can be evaluated analytically because of the simple structure of the two-qubit reduced density 
matrix. A well-known result pertaining to the spin trimer is that there is no pairwise entanglement 
at both T = and at finite temperatures for < e < 1 |18| . We have now shown that the QD has 
a non-zero value for both T = and T ^ 0. An interesting observation is that the QD has a larger 
value, Q (pab) — 0.333, in the FM case than that, Q (pab) — 0.125, in the AFM case. The classical 
correlation, C (pab), has the same value 0.082 in both the cases with magnitude lower than that of the 
QD. On the inclusion of an external magnetic field, the QD jumps in magnitude at the first-order QPT 
point, h c = (figure [3]). This feature is similar to the entanglement jumps seen at first-order QPT 
points [25| 126 j . Dillenschneider |27| has investigated quantum phase transitions in the one-dimensional 
spin- 1 transverse Ising and AFM XXZ models using the QD as a measure. For both the spin models, 
the QD displays behaviour similar to that of entanglement quantified by the concurrence in the vicinity 
of the QPT point. In the case of the AFM spin tetramer, a first order QPT takes place at J x = J 2 
separating the RVB ground states |11), (Ji > J 2 ) and |10), (Ji < J 2 ). Again the QD (as well as the 
classical correlation) exhibits a discontinuity at the transition point. This is true in both the cases of n.n. 
and n.n.n. qubits. A first order QPT involves a discontinuity in the first derivative of the ground state 
energy with respect to a coupling parameter (h in figure [3]) . This implies a discontinuity in one or more 
of the elements of the reduced density matrix at the transition point |26l l2"Tj . Since the QD is dependent 
on the elements of the reduced density matrix, a first order QPT gives rise to a discontinuity in QD. 
In both the cases of the trimer and the tetramer, the asymptotic decay of the QD with temperature 
indicates that thermal fluctuations cannot kill the quantum correlations though the QD is reduced in 
magnitude at higher temperatures. Some recent studies on the thermal QD in spin models [T2j [28] arrive 
at a similar conclusion. 

The successful implementation of quantum computation and communication protocols depends on 
the robustness of entanglement in quantum states. The inevitable interaction between a system and its 
environment results in decoherence and degradation of the entanglement. The entanglement dynamics 
due to decoherence may bring about the complete disappearance of entanglement at a finite time, 
termed the "entanglement sudden death" [29]. Some recent studies [TT f |2" H 130] have shown that the QD, 
in the presence of a Markovian environment (memoryless dynamics), decays in time but vanishes only 
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asymptotically. In fact, Ref. [23] discusses an interesting example of the QD remaining constant up to a 
time t = t with the decay setting in only when t is > t. Our studies of the AFM trimer and the tetramer 
show that the two-qubit reduced density matrices at T = have the form of Werner states. As shown by 
Werlang et al. [12] , the QD vanishes asymptotically with time when an initial Werner state is subjected 
to a dephasing channel. We have demonstrated this for the trimer and the tetramer in Figures ([8]) and 
(J9]) with the 7 — > 1 (7 = 1 — e~ rt ) limit corresponding to t — > 00. The same figures show that the 
pairwise entanglement, as measured by concurrence, is either zero at all times (figure [8]) or undergoes a 
"sudden death'at a finite time (figure [9]). Some recent studies [3T j [32 ] have shown that the use of states 
for which entanglement is zero (mixed separable states) but QD is non-zero, can improve the efficiency 
of certain computational tasks in comparison with classical computing. The spin trimer ground state 
provides an example of a state with zero entanglement and non-zero QD. Ref. [16J provides a number of 
examples of molecular magnets described by spin trimers and tetramers. Molecular spin clusters are ideal 
candidate systems to test quantum information theoretic concepts. Recent advances in supramolecular 
chemistry provide tools to engineer synthetic spin clusters like a molecular cluster of three qubits [33J . 
The role of decoherence causing crossover from the quantum to the classical domain may be ideally 
studied in mesoscopic systems like molecular spin clusters [33] • Since quantum correlations, in terms 
of the QD, persist up to very high temperatures and the decoherence time under specific conditions is 
quite long, such correlations could provide the basis for the implementation of quantum information 
tasks. Molecular magnets, described by small spin clusters, are expected to play an important role in 
such applications. 

In this paper, we have not investigated the dynamics of the QD and the entanglement of the reduced 
two qubit systems subjected to a non-Markovian environment. A number of recent studies [351 1361 1371 138] 
have identified some interesting features of the non-Markovian dynamics of the QD which are absent 
in the Markovian case. The two qubits interact with either independent or common non-Markovian 
environments (reservoirs). In the studies carried out so far, the two-qubits are not coupled to each 
other. For independent reservoirs, the QD is found to vanish only at discrete time points [35} |36| [37] 
whereas the entanglement disappears in a finite time interval. In the case of a common reservoir, the 
entanglement dynamics exhibit damped oscillations whereas the QD is characterised by isolated kinks at 
which there is a jump in its derivatives [37] . Due to the memory effect of the environment, some of the 
quantum correlations lost during the dissipative dynamics can be restored to the qubits giving rise to 
the "sudden birth of entanglemenf'and "revivaf'of the QD [35, 36, [37] [38] . As in the Markovian case [24J, 
for a specific class of initial states, the qubit system exhibits a sudden transition between clssical and 
quantum decoherence at time t = t [38] . For t < t, the amount of classical correlations decays whereas 
the QD remains frozen. The reverse situation holds true for t > t. In the non-Markovian case, multiple 
such transitions can occur due to the memory effect of the environment. For a system of two coupled 
qubits in a non-Markovian environment, there is no study as yet which investigates the dynamics of the 
entanglement and the QD in the same framework. The spin clusters considered in the present paper 
correspond to interacting qubit systems. Studies on the time evolution of quantum correlation in such 
systems interacting with a non-Markovian environment are essential for a fuller understanding of the 
problem of decoherence in molecular magnets. 
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